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ABSTRACT 

Images covering the central 450 x 100" (~ 8.0 x 2.0 kpc) of NGC 5128 (Centaurus A) 
obtained using SCUBA at 850 and 450^m witli beam sizes of 14.5 and 8" respectively, 
are presented. These data are compared with those obtained at other wavelengths, 
in particular the optical, mid-infrared, and far-infrared continuum. The sensitive 850 
and 450/um images show that the submillimeter (submm) continuum morphology and 
spectral index distribution of Centaurus A comprise four regions: an unresolved AGN 
core, an inner jet interacting with gas in the dust lane, an inner disk of radii ~ 90", 
and colder outer dust. The inner disk has a high surface brightness, reverse-S-shaped 
feature in the 850 and 450/im images that coincides with the regions of intense 7 and 
15 //m continuum and a region of active star-formation. The infrared (IR) and submm 
images appear to reveal the same material as predicted by a geometric warped disk 
model consisting of tilted rings. We suggest this scenario is more plausible than that 
recently proposed in literature suggesting that the mid-IR emission in Centaurus A is 
primarily from a bar, with a structure that is different from the extended warped disk 
alone. A dust mass total of 2.2 x lO^M© has been calculated within a radius of 225", 
45% of which is in the star-forming region of radius ~ 90" about the nucleus. 

Subject headings: dust — galaxies: active — galaxies: individual (Centaurus A, NGC 5128) — 
galaxies: structure — submillimeter — radiation mechanisms: thermal 



1. Introduction 

NGC 5128 (Centaurus A), at an assumed 
distance of ~ 3.5 Mpc (Hui et al. 1993), is 
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the nearest giant elliptical galaxy and is re- 
markable in several respects. It is a power- 
ful radio source with a double-lobed structure 
extending approximately 3°. 5 x 8°. 5 on the 
sky (Bolton, Stanley, & Slee 1949; Clarke et 
al. 1992; Tingay et al. 1998). At the other 
extreme of scales a central source, less than 
0.4 milliarcsec (~ 0.008 pc) in extent (Keller- 
mann, Zensus, & Cohen 1997), is prominent 
on radio through X-ray images. This com- 
pact object feeds subparsec-scale relativistic 
outflows that are approximately aligned with, 
and clearly the generators of, the vast outer 
radio lobes (e.g. Tingay et al. 1998). 
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The optical appearance is dominated by 
a dramatic warped dust lane at least 12'. 5 
in E-W extent, which effectively bisects the 
main body of the elliptical galaxy and almost 
completely obscures the nucleus and all op- 
tical structure in the inner 500 parsecs (e.g. 
Schrcier et al. 1996). The outer isophotes of 
Centaurus A are markedly elongated in P.A. 
~ 25°. Faint shells, associated with both 
HI and CO emission, arc evident in these out- 
lying parts (Malin, Quinn, & Graham 1983; Schi- 
monovich et al. 1994; Charmandaris, Combes, & 
van der Hulst 2000). The somewhat chaotic dust 
lane, and especially the shells, are strong evidence 
of a relatively recent merger (Baade & Minkowski 
1954: Graham 1979: Tubbs 1980: Sehreier et al. 
1996; Israel 1998) which is generally believed (c.f. 
Marconi et al. (2000) and references therein) to 
be responsible for the nuclear activity. Ha and 
molecular line observations (e.g. van Gorkom et 
al. 1990) indicate that the nucleus is surrounded 
by a rapidly rotating massive inner disk of ra- 
dius ~2kpc with a pronounced warp (Nicholson, 
Bland-Hawthorne, & Taylor 1992; Quillen et al. 
1992), a scenario that is supported by modelling 
of the structure of the obscuring dust seen in near- 
infrared (Quillen, Graham, & Progel 1993) images. 

Submillimetre wavelength emission, appar- 
ently thermal in origin, was first observed from 
Centaurus A by Cunningham et al. (1986), using 
the single-pixel bolometer UKT14 mounted on the 
3.8-m UK Infrared Telescope with a beam size of 
~ 80". Millimeter- wave continuum observations 
by Eckart et al. (1990) in 22" and 45" beams show 
an unresolved nuclear source, surrounded by ex- 
tended thermal emission seen by IRAS at 50 and 
100/xm and which is roughly co-extensive with a 
region of strong CO molecular line emission. More 
extensive submm continuum observations were 
made by Hawarden et al. (1993), who mapped 
the galaxy using the smaller beamsize of ~ 15" 
afforded by UKT14 mounted on the 15-m James 
Clerk Maxwell Telescope (JCMT).*^ Although ex- 



The James Clerk Maxwell Telescope is operated by 
The Joint Astronomy Centre on behalf of the Parti- 
cle Physics and Astronomy Research Council of the 
United Kingdom, the Netherlands Organization for 
Scientific Research, and the National Research Coun- 
cil of Canada. 



tended emission was visible in their images, these 
were not sensitive enough to reveal much detailed 
structure beyond a general extension correspond- 
ing to the optical dust lane and a brighter, elon- 
gated, central feature, with an apparently thermal 
spectrum, surrounding the strong, flat-spectrum 
nuclear source. 

This paper presents new high-quality images 
of Centaurus A obtained at submm wavelengths, 
that are sensitive in particular to emission from 
cold interstellar dust. These resolve the struc- 
ture of the inner disk down to scales of about 8" 

150 pc), while also revealing details not previ- 
ously seen of the faint outer dust emission at large 
radii. Recently, Mirabel et al. (1999) used ISO- 
CAM, the infrared (IR) imager on the Infrared 
Space Observatory {I SO) J to obtain images at 
mid-IR wavelengths, and these results are com- 
pared, together with observations at other wave- 
lengths, with those presented in this paper. Re- 
processed IRAS^ data are also presented to com- 
plement the detailed study of the dust properties 
in the disk. 

2. Observations 

Centaurus A presents special challenges for 
observations at submm wavelength from Mauna 
Kea since it never rises more than 28° above the 
southern horizon. Images were nevertheless ob- 
tained using the submm wavelength bolometer ar- 
ray SCUBA (Holland et al. 1999) during several 
observing periods with the JCMT during 1998. 

The simultaneous 850 and 450/Ltm images were 
obtained during three nights (UT 1998 March 29, 
April 9 and April 11) during a period in which sky 
conditions were exceptionally transparent, with 
zenith opacities often as low as 0.1 and 0.4 or 
better at 850 and 450/um respectively. The data 
were taken using the "jiggle-mapping" mode of 

^ISO is an ESA project with instruments funded 
by ESA member states (especially the PI countries: 
France (for ISOCAM), Germany, the Netherlands and 
the United Kingdom) with participation of ISAS and 
NASA. 

*The IRAS data were re-processed using HIRES rou- 
tines at NASA/IPAC which is operated by the Jet 
Propulsion Laboratory, Caltech, under contract to the 
National Aeronautics and Space Administration. 
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SCUBA, which provides a dcnscly-oversampled 
image with a spacing of 2". 18, and each observa- 
tion results in a field about 2'. 3 in diameter. The 
array orientation rotates with respect to the sky as 
the observation progresses, further sampling the 
image plane. A series of seven overlapping field 
centers were observed, placed 55" apart along a 
line at a position angle of 120° (North through 
East), roughly corresponding to that of the dust 
lane of Centaurus A. The chop throw was 120" 
perpendicular to this line. 

The strong overlap (more than 50%) between 
neighboring fields compensated for the limited ro- 
tation of Centaurus A with respect to the SCUBA 
arrays, offsetting the effect on the final reduced im- 
ages of occasional excessively noisy pixels. A dif- 
ferent subset of the seven fields was observed on 
each of the three nights. Since the bright nuclear 
source appears on the central three fields, it was 
observed at least twice per night, which allowed 
corrections for pointing drifts to be measured and 
applied to the final images. 

To determine the beam pattern of the JCMT, 
the bright unresolved blazar 3C 279 was observed 
each night using the same jiggle-pattern mode 
as applied to Centaurus A. Atmospheric opac- 
ities were determined from skydips made at in- 
tervals during the observing, and instrumental 
gains were derived from images of the JCMT 
secondary calibrators CRL618, IRC-l-10216 and 
IRAS 16293-2422 (Sandell 1994, 1998). The 
imaging data analysis was undertaken using the 
dedicated SCUBA data reduction software SURF 
(Jenness, Lightfoot, & Holland 1998), as well as 
KAPPA, GAIA and CONVERT software pack- 
ages provided by the Star link Project^. The data 
reduction consisted of first flatfielding the array 
images, and then correcting for atmospheric ex- 
tinction. Next, pixels significantly noisier than the 
mean were blanked-out; and, after initial inspec- 
tion of raw images, pixels containing relatively 
little flux from the source were used to correct for 
correlated sky noise in each individual jiggle-map. 
Corrections for pointing drifts were incorporated 
in the final images using the fields in which the 
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core source appeared. The apparent core source 
size was reduced from 9.24 to 8".42 (FWHM) at 
450/im after the inclusion of these corrections. 

An image at 850/xm had been previously ob- 
tained under somewhat worse conditions on UT 
1998 February 14, 15 and 16. This image was 
made using the "scan-mapping" mode and cov- 
ers an area 4' square centered on the bright core 
source. Despite the lower quality of this dataset, it 
provides a useful comparison dataset for the real- 
ity of features in the two 850/im images. The sky 
transmission was not good enough on this occa- 
sion for the simultaneous 450^m data to be useful. 
The instrumental gains for the "scan-mapping" 
mode observations are 20% less than those typi- 
cal for the "jiggle-mapping" mode; therefore, the 
gains have to be determined separately for the in- 
dividual modes. The beam in the 850/im "scan- 
mapping" image of Centaurus A in this paper is 
slightly larger than in the "jiggle-mapping" im- 
ages and the photometry differences between the 
two final maps are within and proportional to the 
uncertainties: the percentage differences are larger 
in the less sensitive parts of the maps. 

3. Results and Discussion 

3.1. The Submm Images and Compar- 
isons with other Wavelengths 

The images of Centaurus A at 850 and 450/im 
obtained from the combined "jiggle-map" datasets 
are about 6' long and 2' wide (respectively 7kpc 
and l.Tkpc at Centaurus A) and cover much 
of the dust disk. The 850 and 450/xm maps, 
which respectively have rms noise values of 5 and 
20mJy/beam, are shown in Figure 1. Figure 2 
shows the 850 /xm image obtained from the "scan- 
mapping" observations, which has an rms noise of 
25 m.Jy/beam. 

The submm images show extended emission 
centered on the nucleus of Centaurus A and ori- 
ented at a position angle of about 115° (North 
through East) on the plane of the sky, in roughly 
the same orientation as the prominent optical dust 
lane. Within about 90" at both 850 and 450//m, 
the emission has a much higher surface brightness 
than at larger radii. The bright structures have a 
reversed S-shape suggestive of spiral structure or, 
as proposed by Mirabel et al. (1999), a bar. 
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The central source has a flux density per beam 
at least 40 times brighter than the surrounding 
emission from the dust lane, as shown in the pro- 
files along the major axis in Figure 3. Further- 
more, the close resemblance of the central source 
to the profiles of the JCMT beam (line with crosses 
in Figure 3), shows the source is unresolved. The 
point-source core has been isolated from the ex- 
tended emission and separate submm-to-IR spec- 
tral energy distributions (hereafter SED) have 
been constructed for the core and extended galac- 
tic emission. These are discussed in Section 3.3. 

Figure 4 shows the 450/im contours, our highest 
resolution new SCUBA image, superimposed on 
a) an optical waveband (395 and 540 nm) image 
courtesy of the Anglo- Australian Observatory^'^ 
and b) the ISOCAM 7/im image. We now discuss 
these comparisons, as well as what we learn from 
the 450/850^m spectral index distribution map, 
in more detail. 

3.1.1. Submm vs. Optical Morphology 

The well-defined reverse-S-shape structure seen 

in high surface brightness submm emission, within 
90" of the galaxy nucleus, is largely indiscernible 
in optical maps, implying that the denser material 
around the nucleus, seen in the submm- wavelength 
images is heavily obscured at optical wavelengths. 
However, the southern edge of the SE high sur- 
face brightness submm emission aligns with the 
southern ridge of the optical dust lane, in turn sug- 
gesting that the SE reverse-S-shape structure is on 
the near-side and therefore relatively less obscured 
in the optical, as reported by Block & Sauvage 
(2000) and Quillen, Graham, & Frogel (1993) and 
based on their analyses of mid-infrared versus V 
and near-infrared data respectively. Furthermore, 
if the sub mm/optical morphology is a manifesta- 
tion of a warped dusty disk or spiral structure that 
is highly inclined to the plane of the sky (c.f. Sec- 
tion 3.2), the southern edge is consistent with an 
inner fold of the disk seen tangent to the line of 
sight. Such a ridge would have a high column den- 



' Original plate material for the Digitised Sky Sur- 
vey in the southern sky is copyright (c) the Anglo- 
Australian Observatory and was used, with their per- 
mission, to produce the The Digitised Sky Survey at 
the Space Telescope Science Institute under US Gov- 
ernment Grant NAG W-2166. 



sity in the line of sight, as evident in the aligned 
southern edges of high surface brightness submm 
emission and the dark optical lane. 

A very rough estimate can be made of the 
optical depths detected in the SCUBA images, 
allowing a quantitative comparison between the 
submm emission and optical obscuration. Us- 
ing the emissivity (i.e. grain absorption coeffi- 
cient) from Hildebrand (1983) and assuming a 
grain temperature of 17 K (c.f. the diffuse dust 
in the Galaxy), a U-band optical depth of about 
10 can be inferred for a surface brightness of 
~ 30mJy/beam at 850/zm, as seen in the low-level 
emission detected in the current submm images. 
Dust clouds that are evident in the optical image 
presented in this paper are expected to have V- 
band optical depths of ~ 10. Therefore, under 
reasonable assumptions, the dust lanes seen in the 
optical obscuration map exhibited here should be 
detected in the SCUBA images. 

Figure 4 shows that the low surface bright- 
ness submm emission generally follows the optical 
dust absorption distribution, including the clock- 
wise twists in the east and west of the dust lane. 
This submm emission corresponds especially well 
with the dark optical obscuration in the NW and 
northern part of the SE dust lanes, and also with 
the very dark twist in the east, indicating that 
some low level submm emission must arise from 
the optically prominent dust material. There is 
no marked difference between the flux densities of 
submm low level emission that corresponds with 
the dark optical obscuration in the northern edge 
of the SE dust lane and the low level emission in 
the southern ridge of the dust lane: in both re- 
gions, the flux densities are in the range ~ 40 to 
~ 120mJy/beam and ~ 10 to ~ 60mJy/beam at 
450 y^m and 850 /im respectively. Schreier ct al. 
(1996) flnd that the extinction along the southern 
ridge is much more than that in the northern edge 
of the SE dust lane, where it reduces the R flux 
by a factor of 6 compared to a factor of only 1.5 to 
2 in the northern edge of the SE lane. Therefore, 
the similarity of the flux densities in the northern 
and southern low level submm emission despite 
the higher extinction in the southern ridge indi- 
cates that the northern edge of the SE dust lane 
is well in the foreground, in a region of low stellar 
density, and is thus heated by a more dilute stellar 
radiation field. 
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Within the length (major diameter) of 270" 
that excludes the clockwise twists, low surface 
brightness submm emission extends farther north 
and south beyond the optical feature: the average 
width (minor diameter) of the dust feature seen 
in submm emission is >~ 90" compared with the 
optical obscuration of ^< 60" (see Figure 4(a)). 
The extension is appreciably marked in the south, 
where the dust giving rise to the submm emis- 
sion here lies on the far-side of the galaxy and is 
overlain by the stellar body. Consistent with this 
scenario, when Schreier et al. (1996) removed the 
effects of foreground dust obscuration from HST 
I' band images, they found a band of residual ob- 
scuration, presumably caused by dust within the 
stars on the farther side of the galaxy. 

While it is clear that the maximum width at 
which the dust lane will be detected in extinc- 
tion and seen in thermal emission will depend 
on the relative sensitivity of the two techniques, 
the width at which the feature is detected in the 
submm observations (>^ 90", c.f. Figure 4(a)) 
surpasses the maximum width detected in all opti- 
cal observations we arc ciirrcntly aware of, includ- 
ing the recent HST (Sclircicr ct al. 1996; Marconi 
et al. 2000) and VLT" observations (~< 60"), 
supporting the assertion that the low-level submm 
emission extends well beyond the optical obscura- 
tion. 



dust components," attributing the optical dust 
features to "small amounts of cold dust in the 
outer parts" of the system. 

The SCUBA data exhibited in this paper, 
which have much greater sensitivity and extend 
much farther from the nucleus than the images 
published by Mirabel et al. (1999), show that 
there are remarkable similarities between the ap- 
pearance of Centaurus A in the mid-IR and the 
submm (see Figure 1, Figure 4(b)). This is par- 
ticularly true of the 450 and 7/im images, even 
though their wavelengths differ by a factor of 65. 
Both images show the reverse-S-shaped, high sur- 
face brightness striicturc out to 90" from the core 
and the fainter extensions of this structure out 
to 120". The major difference is that the submm 
emission is seen to much larger angular distances 
than is the mid-IR, and as pointed out above, faint 
submm emission from some of the more outlying 
dust in the optical dust lane is evident. 

The contours in Figure 4(b) show a feature ex- 
tending to about 15" radius from the nucleus in 
position angle 145°; this may correspond to that 
attributed to a circumnuclear torus by Hawarden 
et al. (1993); Israel (1998) and Bryant & Hunstead 
(1999). The SCUBA and ISOCAM images suggest 
it may represent the inner folds of the reverse-S- 
shaped structure (or warped disk, c.f. Section 3.2) 
mentioned above. 



3.1.2. Submm vs. ISO Mid-IR Morphology 

Mirabel et al. (1999) presented images of 
Centaurus A at 7 and 15/im obtained with ISO- 
CAM as well as early SCUBA "jiggle- mapping" 
images of a 2' field around the core of the galaxy 
at 450 and 850/im. They note the similarity of 
their IR (ISOCAM) and submm (SCUBA) images 
(which are much less sensitive to faint emission 
than those presented here), observing that there 
was the "same general distribution" between the 
"warm" and "cold" dust seen respectively in the 
ISOCAM and the early SCUBA image published 
by them. Mirabel et al. (1999) conclude that the 
absence of submm emission from the optical dark 
lanes "is not due to major differences between the 
spatial distributions of the cold and very warm 



^^VLT astronomical imago gallery is avail- 
able at http://www.eso.org/outreach/info- 
events /ut Ifl /astroimages.html 



3.1.3. The Spectral Index Distribution 

Using the 450 and 850/im imaging data, the 
global spectral index distribution of Centaurus A 
is derived at submm wavelengths and used to de- 
lineate the submm components in the nuclear re- 
gions of the galaxy as well as to investigate the 
dust properties of the extended emission. The fol- 
lowing procedure was used. The 450/Ltm data were 
first smoothed to the 14" resolution of the 850/im 
map. The spectral index a, where Si, oc i/", be- 
tween 450/im and 850/tm was then computed as 



a = log 



S 



450 



'850 



log 



850 
450 



(1) 



where S is the flux density per beam at each point 
in the map. The uncertainty in the spectral index 
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is then given by 




(2) 

where the variables are as for equation 1. The 
largest probable uncertainty in AS arises from the 
calibration, especially for the 450/Ltm flux densi- 
ties. The maximum uncertainty in our final map 
is estimated as | Aa \^^^^ « 1.5 and the mean un- 
certainty as I Aa lujgan ~ O-^- '^^^ uncertainty in 
the regions with a > 3.2 is most probably under- 
estimated as this is also the region of the map with 
the least sensitivity. 

The spectral index map of Centaurus A be- 
tween 450 and 850/im is shown in Figure 5. There 
are four regions apparent in the map. Contours 
at the spectral indices 2.0 (dashed lines) and 3.0 
and 3.3 (solid lines) are overlayed to guide the eye 
to the less obviously manifested features. In the 
nuclear area, the unresolved core and a feature 
apparently representing the inner jet are distin- 
guished. Farther from the nucleus we can discern, 
with progressively increasing spectral indices, the 
bright, reverse-S-shaped structure, familiar from 
the intensity maps, and the fainter outer dust, 
roughly corresponding to the optical dust lane (see 
Section 3.1.1). 

The spectral index of the unresolved core 
is close to zero, consistent with the nonthermal 
spectrum from an AGN nucleus. This paper con- 
firms that the core source has a flat submm spec- 
trum (i.e. a « 0), which is common for BL Lacs 
(blazars) (e.g. Brown et al. 1989), and as Hawar- 
den et al. (1993) and Kellermann, Zensus, & Co- 
hen (1997) have noted, suggests that Centaurus A 
may harbor a low- luminosity blazar. 

The feature which this paper refers to as the 
inner jet is not evident on the intensity maps 
but is prominent in Figure 5 as (green-blue) areas 
of low spectral index to the NE and SW of the 
nucleus. The spectral index of these features is 
in the range 1.5 < a < 2.5 (see dashed contour 
lines). The northern extension is coincident with 
the inner radio jet (Clarke et al. 1992), which, as 
implied by the designation in this paper, is not 
considered to be merely by chance. However, the 
observed spectral index is much higher than would 
be expected from extrapolation of the power-law 
spectrum of the radio jet itself. One possible ex- 



planation is that the submm spectrum is not a 
manifestation of the electron distribution of the 
relativistic radio jet itself (or the jet overlaid by 
a warm dust component), but rather of some in- 
ternal mechanism (e.g. ionization of gas) resulting 
from the interaction of the jet and the ISM in the 
inner regions of the galaxy. 

Brodie, Konigl, & Bowyer (1983) and others 
have proposed that optical jets in galaxies such 
as Centaurus A may be due to emission from in- 
terstellar gas that has been entrained and heated 
by the flow of relativistic particles from the nu- 
cleus. The entrained gas would generate free- 
free continuTim emission with spectral index, in 
the optically thick case, of a w 2, consistent 
with the observed spectral range of the inner jet 
1.5 < a < 2.5. We attribute the southern nu- 
clear extension in our submm spectral map to the 
same mechanism, though in this case associated 
with the fainter counter jet recently reported by 
Tingay et al. (1998) and Kraft et al. (2000), re- 
spectively from their high resolution VLBI radio 
and Chandra X-ray images. 

A similar interaction of the jet and interstellar 
medium was proposed by Joy et al. (1991) based 
on their near-IR observation of the inner jet of 
Centaurus A, which manifested itself as a region 
of bluer color, with spectral index 1.3 < a < 2.9, 
coincident with the radio jet. Marconi etal. (2000) 
also reported a 'blue channel' coincident with the 
radio jet, but they attributed it to low extinction 
resulting from relatively low concentration of dust 
that has been mechanically 'evacuated' by the jet. 
This paper disagrees with this interpretation as 
there is no sign of a deficiency of emission on the 
450/Ltm map at this location; the inner jet feature 
of Figure 5 arises from a relative enhancement of 
emission at 850/im in these locations. 

The spectral index of the inner disk is in the 
range 2.5 < a < 3.2, consistent with thermal emis- 
sion from dust. The lowest indices occiir in an al- 
most circular region about 30" in radius around 
the nucleus; this is coincident with an area of hot 
molecular gas reported by Schreier et al. (1998)^^ 
from their HST observations. Emission with a 
steeper spectrum (2.6 < a < 3.2, see solid con- 



STScI-PRC98-14 Electronic Press Release is available 
at http://oposite.stsci.edu/pubinfo/pr/1998/14/pr- 
photos.html 
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tour lines) approximately traces the inner disk of 
brighter structures seen in the 7 and 450/Ltm im- 
ages (see Section 3.1.2 and Figure 4b) out to ~ 
90". Farther still from the nucleus, in the fainter 
outer disk the spectral index increases to a 4. 
As underscored by the solid line contours in Fig- 
ure 5, the spectral index is asymmetric about the 
nucleus in the inner disk region of the map: it is 
somewhat flatter on the eastern part of the inner 
disk, indicating slightly different relative distribu- 
tions of dust grains on either side of the nucleus. 

If the dust seen in the maps presented in 
this paper is emitting with spectral index a = 
2 + 13 (where [3 is the cmissivity index) , the spec- 
tral index in the extended region implies that for 
a given temperature, the dust in Centaurus A has 
on average a low /?., i.e. the dust is made of rel- 
atively large grains, and that its temperature is 
fairly cool. Furthermore, the somewhat flatter 
spectral index in the south-east suggests the dust 
in this region is made of smaller grains with mod- 
erately warmer temperatures than those in the 
same region in the NW. This is not surprising as 
the SE is reported to exhibit stronger H I emission 
(van Gorkom et al. 1990), weaker ^^00(1-0) and 
stronger Ha line emission compared to the NW, 
that led Eckart et al. (1990) to speculate that star- 
formation and its associated high radiation field 
were greater in the SE. 

The uncertainty in the emissivity analysis in 
our spectral index map (Figure 5) does not per- 
mit us to quantitatively examine the effects, if any, 
due to diffuse "cirrus" grains in the map. The 
dust in the outer disk, where the emissivity in- 
dex j3 > 1.5, may well consist of "cirrus" grains, 
which, as it is generally accepted (e.g. Rowan- 
Robinson 1992), would have an emissivity index 
close to /3 = 2 and temperatures 15 < T < 40 K. 
Eckart et al. (1990) associated the far-IR emis- 
sion outside the molecular star-forming disk (ie. 
radius >^ 90") with "cirrus" clouds with scale 
height larger than that of the molecular gas disk. 
They predicted that the submm extended emission 
would also originate from dust constituting "cir- 
rus" grains as could well be the case. The large 
western and eastern warps seen in the 850jum and 
450/im images (see Figure 1) resemble the distri- 
bution of the H I emission (Schimonovich et al. 
1994), indicating that, if the dust in the warp fea- 
tures (and the general outer disk) in Centaurus A 



is composed of "cirrus" grains, the H I emission 
traces the "cirrus" dust in NGC5128 just as in 
the Milky Way. 

3.2. Implications of the Multiwave- 
length Images 

Our submm images demonstrate that outly- 
ing dust in Centaurus A is significantly cooler 
than the material in the bright elongated features 
within ~ 90" of the nucleus. This is naturally 
to be expected: at larger radii the ambient radia- 
tion field heating the dust is more dilute because 
of the much lower density of stars (Eckart et al. 
1990; Schreier et al. 1996; Marconi et al. 2000). 
The emission from this dust generally follows the 
optical dust lane. 

The warped disk model proposed by (Quillen 

et al. 1992) and explored by Quillen, Graham, & 
Frogel (1993), consisting of tilted rings of material, 
predicts the structure of such warmer and colder 
material, respectively seen in the mid-IR and far- 
IR to submm, rather well. This scenario is sup- 
ported by Eckart, Wild, & Ageorges (1999), who 
showed that a warped structure of tilted rings ex- 
plained not only the line emission in the disk of 
Centaurus A (Quillen, Graham, & Frogel 1993), 
but also a complex absorption- line system towards 
the nucleus, without a need for any additional 
structures. 

From their ISO 7 and 15/im images, Mirabel 
et al. (1999) conclude from a comparison with 
similar mid-IR images of the dwarf barred spi- 
ral NGC 1530 that the structure imaged in the 
mid-IR in Centaurus A is itself a barred spiral: 
the prominent emission peaks ~ 75" from the nu- 
cleus at each end of the bright mid-IR and submm 
structures are hypothesised to be foreshortened 
arms twisting anti-clockwise from the outer ends 
of the bar. Block & Sauvage (2000) support this 
proposal from their VHK and 15/im study of 
Centaurus A, noting too that the warm dust seen 
by ISOCAM contributes little to the extinction 
seen at optical wavelengths. 

Kinematics is clearly important in understand- 
ing the true nature of the near-nuclear structures 
in Centaurus A. Mirabel ct al. (1999) address 
these in their Figure 4. a comparison of the IR 
images with the CO kinematics in the form of 
a position-velocity (P-V) plot from Quillen et al. 
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(1992). On this figure they identify the strong CO 
feature extending from —1 .2, 300 kms to +1 .2, 
800kms~^ as the bar, undergoing sohd body rota- 
tion, and weaker outfying features at nearfy con- 
stant velocity at larger radii as the "typical" rota- 
tion curve of galactic disks. 

However, an alternative interpretation of the 
PV diagram is that the "bar" represents a nearly- 
complete, highly foreshortened, ring of material 
at a radius of ~ 65" (~ 1300 pc). In that scenario 
the concentrations of IR and submm emission NW 
and SE of the nucleus represent the ends of the 
ring where the optical depth is maximized. Sim- 
ilarly, the "EW high-velocity feature" mentioned 
by Mirabel et al. (1999), which they note does not 
fit their scenario, is readily explained by another 
ring of molecular material at a radius of ^ 20" 
(~ 400 pc), tilted relative to the outer ring and 
evident on the mid-IR images as two small exten- 
sions to the nuclear feature in PA 80° and 260°. 
Furthermore, the outer features of the PV diagram 
also have an alternative interpretation: the fea- 
ture >1 .7SE of the nucleus has the properties of 
an arc of material starting at radius of ~ 4.6 kpc 
and ending perhaps 2 kpc from the nucleus close to 
the minor axis: in an normal edge-on spiral, such 
a feature would be interpreted as a spiral arm. 

Marconi et al. (2000) also accept the bar model 
since its edges appear to be delineated by linear 

concentrations of star- forming regions seen in their 
Pa a images, suggesting that these delineate the 
shocks normally seen as dust lanes along the lead- 
ing edges or the midlines of galaxy bars in optical 
images. However, they also noted that their data 
do not rule out a warped ring without a bar, as the 
star-formation is found in regions of the warp that 
are tangent to the line of sight, and if the young 
stars are above or inside the disk, they would in- 
deed appear to be concentrated where the Pa a 
emission is seen. 

We suggest the latter interpretation is more 

plausible for the following two reasons. First, the 
same elongated distribution of Pa a emission re- 
gions, which Marconi et al. (2000) show to coin- 
cide precisely with the elongated mid-IR ridge that 
passes just south of the nucleus, is equally coinci- 
dent with the dense, narrow, optical dust lane that 
passes south of the nucleus in Figure 4(a). Since 
this is the only feature of the bright central com- 



plex (apart from the nucleus itself) that is poten- 
tially identifiable with an optically-visible feature, 
we suggest that it is not in fact part of the near- 
nuclear complex at all, but, as suggested by the 
models of Quillen et al. (1992), a manifestation of 
an outlying fold or ring in the warped disk. 

Second, we note that the bar-shocks traced 
by dust in early-type field spirals are asymmet- 
ric about their nuclei, since they lie in the leading 
edges of the bars, while those in late type systems 
arc centered in the bars. However, the Pa a emis- 
sion regions in Figure 13 of Marconi et al. (2000) 
pass by several arcsecs south of the nucleus from 
the SE to west of the nucleus. This is not what 
would be expected of a shock in a bar in a field 
spiral, but more the distribution to be expected 
if the star-forming regions are associated with the 
trailing dust lanes in a typical spiral arm, or in 
this case the dusty rings postulated by Quillen et 
al. (1992). 

Therefore, the combined evidence to date from 
the CO kinematics, the mid-IR images, the bright 
and faint submm features and the Pa a images 
supports the warped rings model of Quillen et al. 
(1992) rather than indicating the presence of a 
true bar, and may even represent spiral structure 
in the dust lane of Centaurus A. 

One of the attractions of the bar scenario for 
both Mirabel et al. (1999) and Marconi et al. 
(2000) is its expected utility as a mover of mate- 
rial from larger to smaller radii in the disk in order 
to fuel the AGN. However, since the association of 
barred structure with the presence of an active 
nucleus in field galaxies, though long sought, is 
marginal at best (c.f., for example. Ho, Fillipenko, 
& Sargent (1996) and references therein), the pres- 
ence of the AGN cannot be taken as an argument 
for the reality of the putative bar. Not having a 
bar to fuel the active nucleus is not too great a 
loss; e.g. recent work by Duschl, Strittmatter, & 
Biermann (2000) suggests that viscosity in a thin 
gas plane may also provide an efficient AGN fu- 
elling mechanism. 

3.3. The Extended Emission Temper- 
ature and Dust Mass Estimates 

This paper has shown that the central source in 
Centaurus A is unresolved and has a flux density 
per beam at least 40 times brighter and a spec- 
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trum markedly flatter than the surrounding emis- 
sion from the dust lane; i.e. the central source is 
clearly distinct from the extended emission (c.f. 
Sections 3 & 3.1.3 and Figure 3 & Figure 5). 
Therefore, an analysis of the SED for the ex- 
tended submm to infrared emission in this paper 
is preceded by the isolation of the point-source 
core from the extended emission and the deter- 
mination of the respective flux density estimates 
for these separate components. Adjustments are 
made when dealing with the JCMT and ISO 
data, both taken with small beams, and the IRAS 
data, which were obtained with significantly larger 
beams. 

The core flux densities from 850 to 7/im are pre- 
sented in Table 1, together with integrated fiux 
densities determined for two regions of the ex- 
tended emission from the dust lane: (1) an ellipse 
of 60 X 180" minus the core measurement and (2) 
an elliptical annulus of 120 - 60 x 450 - 180". The 
far-IR flux densities are archival IRAS data for 
Centaurus A that we re-analysed using the HIRES 
facility to extract the best possible spatial resolu- 
tion and to ensure that we have the same registra- 
tion of the regions from which we determine inte- 
grated flux densities. The ellipse and annulus are 
denoted respectively as the inner disk and outer 
disk, as the two regions have different submm con- 
tinuum morphology and spectral index distribu- 
tions (see Sections 3.1.2 & 3.1.3). The SED for the 
two extended regions are plotted in Figure 6, and 
temperatures as well as dust masses determined 
for the emission from these regions are discussed 
below. 

The compact core is unresolved on the submm 
and ISO mid-IR images (FWHM < 8" at 450/im, 
< 4" at 7/um), as would be expected if it is pre- 
dominantly a manifestation of the unresolved cen- 
tral radio source (c.f. Kellermann, Zensus, & Co- 
hen 1997)). In the case of the IRAS data, the 
best HIRES beams achieved in this paper are still 
fairly large (50" to 90") so that the core extent 
is impossible to constrain and its flux density not 
easily separable from that of the extended emis- 
sion. Therefore, the ellipse IRAS flux densities 
listed in the table are simply the integrated mea- 
surements minus the core flux density of 8 ± 2 Jy, 
at both 100/xm and 60/xm, which is an extrapola- 
tion from the submm core SED (affc^^^ = -0.01, 
S^ocv") to 60/im (5000 GHz). A full discussion of 



the SED of the submm core in Centaurus A, and 
its implications, is outside the scope of this pa- 
per; for an extensive analysis of this matter, sec, 
among many others. Bailey et al. (1986); Hawar- 
den et al. (1993); Kellermann, Zensus, & Cohen 
(1997); Alexander et al. (1999); and Marconi et 
al. (2000). 

The integrated submm to far-IR flux densities 
are flt by a two-component, optically-thin modi- 
fied blackbody of the form 

5, = (1 - exp{-XJXf)[nMTi) + 02-B.(T2)], 

(3) 

where Si, is the observed flux density at frequency 
v, fl the solid angle for the modified blackbody 
component, B^{T) the Planck function at tem- 
perature T, Ao the wavelength at which the op- 
tical depth is unity and /3 the emissivity index of 
the grains. First, a 12K and 40 K two-component 
modified blackbody with dust emissivity index of 
1.3 fits the data for the inner disk, and secondly, a 
12 K and 30 K modified blackbody with dust emis- 
sivity index of 1.6 fits the data for the outer disk 
(c.f. Figure 6). These fits, in particular the de- 
rived relative emissivities and temperatures, are 
consistent with the spectral index maps, showing 
that the average dust temperature decreases from 
the galactic core outwards. Unsurprisingly, these 
temperatures are also consistent with the 42 K and 
32 K determined by Eckart et al. (1990) from their 
50 and 100/im IRAS flux densities alone, as well 
as the average ~ 40 K determined by Unger et al. 
(2000) from their 40 to lOO/xm ISO LWS maps of 
the inner disk region. 

The mass of emitting dust can then be es- 
timated from (e.g. Hildebrand 1983) 

where Si, is the measured flux density at frequency 
i^, D is the distance to the source, B{y,T) the 
Planck function and k^ = ZQ^/'^ap the grain mass 
absorption coefficient where a and p are respec- 
tively the grain radius and density. Values of 
^450 /xm ^ o.25m2kg-i and A:^°'^'" = S.Sm^kg-^ 
(e.g. Hildebrand 1983) are assumed, yielding dust 
masses of 2.5 x lO^M© for T = 40K and 9.6 x 
lO^M© for T = 12K in the inner disk and dust 
masses of 4.9 x lO^Mg for T = 30K and 1.2 x 
lO^Mn for T = 12K in the outer disk. The to- 



9 



tal mass for the dust that emits from the far- 
IR through the submm wavelength is then 9.9 x 
1O^M0 and 1.2 x lO'^M©, respectively in the in- 
ner disk and outer disk, giving an overall total of 
2.2 X IO^Mq for the dust lane region of the galaxy 
we have observed (ie. radii < 225" or ~ 4.5 kpc). 

The advantage of using optically thin submm 
emission to determine the dust mass is, \mfortu- 
nately, offset by the increased uncertainty in the 
properties of dust and subsequently of /cd, as A is 
increased from the far-IR to submm wavelengths. 
A different choice of fed could result in an estima- 
tion of the dust mass that differs by a factor as 
large as ^ 10 (Draine 1990). The dust mass de- 
rived in this paper in a region of radius ^ 90" (the 
inner disk) is at the lower end of the dust mass 
of 1 - 2 X IG^Mq derived by Block & Sauvage 
(2000) from their V — 15y^m spatial distribution 
map in the same area. This is surprising because 
they reported the estimate was from dust in the 
"bar-|-arms of a mini-spiral, excluding any diffuse 
dust." However, given the uncertainty, the mass 
in this paper for the inner disk dust is at least 
consistent with Block & Sauvage (2000). 

4. Conclusions 

The SCUBA 850 and 450/um images in this paper 
show that the submm continuum morphology and 
spectral index distribution of Centaurus A com- 
prise: 

1. The Nucleus and associated structures: 

A distinct, unresolved, flat-spectrum AGN 
core with circumnuclear structure, includ- 
ing, to the northeast and southwest, areas 
of low spectral index which we suggest arise 
from free-free emission in ionised gas en- 
trained in the nuclear outflow; 

2. Inner disk: A prominent, elongated fea- 
ture that may be a circumnuclear ring (or 
bar) extends across the center of the galaxy 
out to a radius of ~ 90". It contains con- 
siderable real structure, including concen- 
trations of emission towards its ends and 
a reverse-S-shaped twist; this structure ap- 
pears reflected about the nucleus. The de- 
tails of this structure are very similar indeed 
to what is seen on mid-IR ISO images, so 
the cooler dust seen in the sub-mm is very 



closely co-extensive with the much warmer 
material (transiently heated by UV-photon 
absorption) seen at 7 and 15 /um, and which 
establishes that the inner disk is a locus of 
vigorous star formation. It is also a strong 
source of CO emission with kinematics con- 
sistent with a ring or rings of gas around 
the nucleus, rather than a bar, as has been 
suggested elsewhere (Mirabel et al. 1999). 

3. Outer disk and outlying dust features: 

Low-k;vel emission, at least some of which 
comes from the foreground material in the 
prominent optical dust lane, extends to the 
eastern, and probably to the western, end 
of our maps. It traces the clockwise twist of 
the optical feature. 

Some of the observed low-level emission 
does not coincide with optically-visible dust. 
The sensitive SCUBA images are, for the 
first time, showing the direct detection of 
emission arising from dust that is in the far 
side of the galaxy, overlaid and obscured 
by the stellar component of NGC 5128, and 
thus not seen in optical images. 

A warped disk model consisting of tilted rings 
(Quillen, Graham, & Frogel 1993) predicts the 
structure of the warmer and colder material, es- 
pecially in the inner disk, rather well. It appears 
that in this vicinity the IR and submm images 
to a large extent reveal the same material; how- 
ever, at larger radii, the dust is cooler because it 
is immersed in a much less intense stellar radi- 
ation field. Alternative arguments are presented 
in papers by Mirabel et al. (1999) and Block & 
Sauvage (2000), who interpret the mid-lR (and 
brighter submm) structures as a true barred spi- 
ral. 

Using the continuum integrated flux densi- 
ties from far-IR through submm wavelengths, 
this paper derives a total dust mass of 2.2 x 
10®M,5 within a radius of 225" of the nucleus 
of Centaurus A. About 45% of the dust mass is 
in the star-forming inner disk within about < 90" 
of the nucleus. 
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Fig. 1. — Images of Centaurus A obtained 
with SCUBA at 850;um (top) and 450/im 
(bottom), derived from a series of individual 
"jiggle- mapping" (see text) exposures offset 
along a line roughly corresponding to the op- 
tical dust band of Centaurus A. The 450/im 
image is smoothed slightly and the central 
source displayed saturated in order to high- 
light the low-level, extended emission. The 
panel key indices are in Jy/beam correspond- 
ing to color-coded intensities, and the contour 
heights in the 850 and 450/Ltm images are [0.01, 
0.03, 0.05, 0.07, 0.1, 0.14, 0.2, 0.3, 0.4, 1.0, 
2.5, and 4.0 Jy/beam] and [0.04, 0.16, 0.28, 
0.40, 0.52, 0.8,1.0, 1.5, 2.5, and 3.5 Jy/beam] 
respectively. 

Fig. 2. — Image of the central 4' at 850//m 
obtained using with SCUBA and using the 
complementary "scan-map" technique. The 
central source is displayed saturated in or- 
der to highlight the low-level, extended emis- 
sion. This map covers a smaller area than 
the "jiggle- mapping" set in Figure 1. The 

This 2-column preprint was prepared with the AAS 
macros v5.0. 



panel key indices are in Jy/beam correspond- 
ing to color-coded intensities, and the contour 
heights are [0.07, 0.1, 0.2, 0.4, 1.0, 2.5, and 
4.0 Jy/beam). 

Fig. 3. — Profiles of the emission from 
Centaurus A at 850 (top) and 450/um (bot- 
tom) along an axis roughly coincident with 
the optical dust lane. The profile of the JCMT 
beam is superimposed in a line marked with 
crosses. This beam profile is obtained from 
a map of the JCMT pointing source 3C 279 
scaled to the peak flux of CentarusA and 
determined along an axis roughly coincident 
with the major axis of the dust lane. 

Fig. 4. — Contours of the 450/xm image 
(see Figure 1) superposed on, from top to 
bottom, (a) an optical waveband image of 
Centaurus A, digitised from the IllaJ emul- 
sion photographic plates which are sensitive 
to emission between 395 and 540 nm wave- 
lengths, courtesy of the Anglo- Australian Ob- 
servatory, and (b) a 7fim ISOCAM negative 
image (Mirabel et al. 1999). The contour 
heights of the 450/im emission in (a) and (b) 
are [0.04, 0.16, 0.28, 0.40, 0.52, 0.8, 1.0, 1.5, 
2.5, and 3.5 Jy/beam.] and [0.4, 0.7, 0.96, 1.2, 
2.4, and 4.8 Jy/beam], respectively. In (a), 
the contours are of the slightly smoothed im- 
age, presented as the 450/xm map in Figure 1 
to highlight the low surface brightness emis- 
sion that corresponds with the outlying parts 
of extended optical dust lane. 

Fig. 5. — Map of the spectral index distribu- 
tion at submm wavelengths, derived from the 
450 and 850/xm data presented in this paper 

(sec Figure 1). Four regions are clearly ap- 
parent in the map: the unresolved core (blue- 
black) and a feature apparently representing 
the inner jet (green-blue) in the nuclear area; 
and the bright, reverse-S-shapcd inner disk 
(yellow-grccn-rcd), familiar from the intensity 
maps (shown in Figure 1), and the fainter 
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outer disk (rcd-ycUow) farther from the nu- 
cleus. Contour heights at the spectral indices 
2.0 (dashed lines) and 3.0 and 3.3 (solid lines) 
are overlayed to guide the eye to the the in- 
ner jet and inner disk features in the map. 
The panel key indices correspond to color- 
coded spectral indices obtained between 850 
and 450/Ltm. 

Fig. 6. — Integrated flux densities in two an- 
nuli centered on the core of Centaur us A. The 
error bars indicate the uncertainty in the flux 
density measurements, which are dominated 
by the uncertainty in the flux calibration and 
to a lesser extent the beam error-lobe contri- 
bution. The top and bottom SEDs are respec- 
tively for the inner disk of 60 x 180" minus 
the core flux density and the outer disk annu- 
lus of 120 - 60 X 450 - 180". The submm to 
far-IR data points arc fit by two component 
optically thin modified blackbodies of temper- 
atures 12 K and 40 K with dust emissivity in- 
dex of 1.3 for the inner disk (top) and 12 K 
and 30 K with dust emissivity index of 1.6 for 
the outer disk annulus (bottom). 
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Table 1: Submm to Mid-IR Fluxes of the 
Centaurus A Core and Extended Emission 



Filter 


Frequency 


Core^ 


Inner disU^ 


Outer dislf 




(GHz) 


(Jy) 


(Jy) 


(Jy) 


SCUBA 850^ 


350 


8.1 ± 0.8 


2.7 ± 0.4 


3.4 ± 0.5 


SCUBA 750^^ 


407 


8.1 ± 1.6 






SCUBA 450^1 


667 


7.9 ± 0.8 


16.9 ± 3.2 


22.6 ± 5.2 


SCUBA 350^ 


866 


7.7 ± 1.9 






IRAS 100 


3000 




119 ± 20*^ 


181 ± 36 


IRAS 60 


5000 




96 lb 15'^ 


77 ± 15 


IRAS 25 


12000 




15 ± 2^ 


6 ± 1 


IS0CAM15 


20000 


1.2 


10.4 




IRAS 12 


25000 




11 ± 2^ 


6 ± 1 


IS0CAM7 


42857 


0.7 


9.4 





"Peak fluxes of the unresolved core, which in the case of 
the SCUBA and ISOCAM images was clearly distinct 
from the extended emission. 

'The integrated fluxes determined in an ellipse of 60 x 
180" minus the core flux. 
■^The integrated fluxes determined in an elliptical an- 

nulus of 120 - 60 X 450 - 180". 

"^The SCUBA filter bandwidths are all about 30 GHz. 

•^In the case of these 100 and 60/Ltm ellipse fluxes, the 
core fluxes subtracted off are an extrapolation of the 
submm core flat SED {af^^^t^ = -0.01, S^ oc v") to 
60 /im (5000 GHz), as the HIRES beams are very large 
(50" to 90") so that it was impossible to isolate the 
core measurements from the extended emission. 

■'^Corc fluxes have not been subtracted from these in- 
tegrated fluxes as no reasonable estimate of the core 
measurements at these wavelengths could be extrapo- 
lated from our present data. 

Note. — The IRAS data were re-processed using HIRES routines at NASA/IPAC. 
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